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Abstract 

We consider a quantum system S interacting sequentially with independent systems 
£rm m — 1,2, .. . Before interacting, each £m is in a possibly random state, and each 
interaction is characterized by an interaction time and an interaction operator, both 
possibly random. We prove that any initial state converges to an asymptotic state 
almost surely in the ergodic mean, provided the couplings satisfy a mild effectiveness 
condition. We analyze the macroscopic properties of the asymptotic state and show 
that it satisfies a second law of thermodynamics. 

We solve exactly a model in which S and all the Sm are spins: we find the exact 
asymptotic state, in case the interaction time, the temperature, and the excitation 
energies of the Sm vary randomly. We analyze a model in which 5 is a spin and the 
£m are thermal fermion baths and obtain the asymptotic state by rigorous perturbation 
theory, for random interaction times varying slightly around a fixed mean, and for small 
values of a coupling constant. 

1 Introduction 

This paper is a contribution to rigorous non-equilibrium quantum statistical mechanics, 
examining the asymptotic properties of random repeated interaction systems. The paradigm 
of a repeated interaction system is a cavity containing the quantized electromagnetic field, 
through which an atom beam is shot in such a way that only a single atom is present in 
the cavity at all times. Such systems are fundamental in the experimental and theoretical 
investigation of basic processes of interaction between matter and radiation, and they are 
of practical importance in quantum optics and quantum state engineering [15, 16, 17]. 

A repeated interaction system is described by a "small" quantum system S (cavity) 
interacting successively with independent quantum systems £i, £2,- ■ ■ (atoms). At each 
moment in time, S interacts precisely with one Sm (with increasing index as time increases), 
while the other elements in the chain C = £1+^2+' ■ ■ evolve freely according to their intrinsic 
(uncoupled) dynamics. The complete evolution is described by the intrinsic dynamics of S 
and of £m, plus an interaction between S and Ern, for each m. The latter consists of an 



*CNRS-UMR 8088 and Departement de Mathematiques, Universite de Cergy-Pontoise, Site Saint-Martin, 
BP 222, 95302 Cergy-Pontoise, France. Email: laurent.bruneau@u-cergy.fr, http://www.u-cergy.fr/bruneau 

tinstitut Fourier, UMR 5582, CNRS-Universite de Grenoble I BP 74, 38402 Saint-Martin d'Heres, France. 
Partly supported by the Mimstere frangais des affaires etrangeres through a sejour scientifique haut niveau; 
Email: Alain. JoyeOujf-grenoble.fr 

^Department of Mathematics and Statistics, Memorial University, St. John's, NL, AlC 5S7, Canada. 
Partly supported by the Mimstere frangais des affaires etrangeres through a sejour scientifique haut niveau; 
Email: merkli@math.mun.ca, http://www.math.mun.ca/~merkli/ 



interaction time Tm > 0, and an interaction operator Vm (acting on S and Sm)', during the 
time interval [ti + • • • + r^-i, ti + • • • + Tm), S is coupled to £m via a coupling operator 
V^. One may view C as a "large system", and hence 5 as an open quantum system. 
Prom this perspective, the main interest is the effect of the coupling on the system S. 
Does the system approach a timc-asymptotic state? If so, at what rate, and what are the 
macroscopic (thermodynamic) properties of the asymptotic state? Idealized models with 
constant repeated interaction, where £jn = £, = T,Vjn = V, have been analyzed in [7, 17]. 
It is shown in [7] that the coupling drives the system to a r-periodic asymptotic state, at an 
exponential rate. The asymptotic state satisfies the second law of thermodynamics: energy 
changes are proportional to entropy changes, with ratio equal to the temperature of the 
chain C. In experiments, where repeated interaction systems can be realized as "One- Atom 
Masers" [15, 16, 17], S represents one or several modes of the quantized electromagnetic field 
in a cavity, and the £ describe atoms injected into the cavity, one by one, interacting with 
the radiation while passing through the cavity, and then exiting. It is clear that neither 
the interaction (r^, Vm), nor the state of the incoming elements £m can be considered 
exactly the same in each interaction step m. Indeed, in experiments, the atoms are ejected 
from an atom oven, then cooled down before entering the cavity - a process that cannot be 
controlled entirely. It is therefore natural to build a certain randomness into the description. 
For instance, we may consider the temperature of the incoming £ or the interaction time r 
to be random. (Other parameters may vary randomly as well.) We develop in this work a 
theory that allows us to treat repeated interaction processes with time-dependent (piecewise 
constant) interactions, and in particular, with random, interactions. Wc are not aware of 
any theoretical work dealing with variable or random interactions, other than [8]. Moreover, 
to our knowledge, this is the only work, next to [8], where random positive temperature 
Hamiltonians (random Liouville operators) are examined. 
The purpose of the present paper is twofold: 

- Firstly, we establish a general framework for random repeated interaction systems and 
we prove convergence results for the dynamics. The dynamical process splits into a decaying 
and a fiucutating part, the latter converging to an explicitly identified limit in the crgodic 
mean. To prove the main convergence result, Theorem 1.2 (see also Theorems 3.2 and 3.3), 
we combine techniques of non-equilibrium quantum statistical mechanics developed in [7] 
with techniques of [8], developed to analyze infinite products of random operators. We 
generalize results of [8] to time-dependent, "instantaneous" observables. This is necessary 
in order to be able to extract physically relevant information about the final state, such as 
energy- and entropy variations. We examine the macroscopic properties of the asymptotic 
state and show in Theorem 1.4 that it satisfies a second law of thermodynamics. This law 
is universal in the sense that it does not depend on the particular features of the repeated 
interaction system, and it holds regardless of the initial state of the system. 

- Secondly, we apply the general results to concrete models where 5 is a spin and the £ 
are either spins as well, or they are thermal fermion fields. We solve the spin-spin system 
exactly: Theorem 1.5 gives the explicit form of the final state in case the interaction time, 
the excitation level of spins £ or the temperatures of the £ are random. The spin-fermion 
system is not exactly solvable. We show in Theorem 7.1 that, for small coupling, and 
for random interaction times r and random temperatures /? of the thermal fermi fields £, 
the system approaches a deterministic limit state. We give in Theorem 1.6 the explicit, 
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rigorous expansion of the limit state for small fluctuations of r around a given value tq. 
This part of our work is based on a careful execution of rigorous perturbation theory of 
certain non-normal "reduced dynamics operators" , in which random parameters as well as 
other, deterministic interaction parameters must be controlled simultaneously. 

1.1 Setup 

The purpose of this section is to explain parts of the formalism, with the aim to make our 
main results, presented in the next section, easily understandable. 

We first present the deterministic description. According to the fundamental principles 
of quantum mechanics, states of the systems S and 6m are given by normalized vectors (or 
density matrices) on Hilbcrt spaces TCg and TCs^, respectively. We assume that dimTCs < 
oo, while the Tis^ may be infinite dimensional. Observables of S and £m are bounded 
operators forming von Neumann algebras dJts C BiJ-Ls) and '^Cm ^ ^0~^£m)- Observables 
As G and As^ G evolve according to the Heisenberg dynamics M. 3 t ag{As) 

and M 9 t a\^{A£^) respectively, where and a^-^ are *-automorphism groups of 9K5 
and OJtf^, respectively, sec e.g. [5]. The Hilbert space of the total system is the tensor 
product Ti = TLs ® "^Ci where Tic = ®m>i '^£m is the Hilbert space of the chain, and the 
non-interacting dynamics is defined on the algebra SDt^ (S)m>i by Oig (S)m>i ■ "^^^ 
infinite tensor product Ti. is taken with respect to distinguished "reference states" of the 
systems S and £m, represented by vectors tps G Tts and tjj^^^^ G H^-,,/. Typically, one takes 
the reference states to be equilibrium (KMS) states for the dynamics a^, a^^, at inverse 
temperatures Ps, P£m- 

It is useful to consider the dynamics in the Schrddinger picture. For this, we implement 
the dynamics via unitaries, generated by self-adjoint operators Lg and L^^, acting on 
BiTis) and BiTLs^), respectively. The generators, called Liouville operators, are uniquely 
determined by 

=e"^#A#e-'*-^#, tGR, and L^V'* = 0, (1.1) 

where # stands for either S or £m ^. In particular, (1.1) holds if the reference states are 
equilibrium states. Let > and V,n G '^s®'^£rn the interaction time and interaction 
operator associated to S and £m- We define the (discrete) repeated interaction Schrodinger 
dynamics of a state vector ip e H, for m > 0, by 

f/(m)^ = e"'^'"-^"' • • • e-'^2L2g-iriLi^^ ^^ 2) 

where 

Lk = Lk + ^L£^ (1.3) 

^ Those vectors are to be taken cyclic and separating for the algebras ftHs and OTg^, respectively [5]. 
Their purpose is to fix macroscopic properties of the system. However, since dimHs < 00, the vector ips 
does not play any significant role. In practice, it is chosen so that it makes computations as simple as 
possible. 

^The existence and uniqueness of satisfying (1.1) is well known under general assumptions on the 
reference states V'# [5]. 
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describes the dynamics of the system during the time interval [ri + • • • + Tk-i, ti + • • • + r^), 
which corresponds to the time step k of the discrete process, with 

Lk = Ls + L£^ + Vk, (1.4) 

acting on Hs (8> 'He,,. (We understand that the operator L^^ in (1.3) acts nontrivially only 
on the n-th factor of the chain Hilbert space He-) 

An operator p on TC which is self-adjoint, non-negative, and has unit trace is called a 
density matrix. A state q{-) = Tr(p-), where Tr is the trace over 7i, is called a normal 
state. Our goal is to understand the large-time asymptotics (m oo) of expectations 

g {U{myOU{m)) = ^(a'"(0)), (1.5) 

for normal states g and certain observables O. Important physical observables are rep- 
resented by operators that act either just on S or ones that describe exchange processes 
between S and the chain C. The latter are represented by time-dependent operators because 
they act on S and, at step m, on the element £m which is in contact with S. We define 
instantaneous observables to be those of the form 

= As'^'j=_iB(j,\ (1.6) 

where G SJt^ and Bm^ G dJt^^^. (we do not write identity operators in the tensor 
product). The class of instantaneous observables allows us to study all properties of S 
alone, as well as exchange properties between S and C. 

Let us illustrate our strategy to analyze (1.5) for the initial state determined by the 
vector ipQ = tjjs ^ ipc, where ipc = ^m>ifp£m- We use ideas stemming from the algebraic 
approach to quantum dynamical systems far from equilibrium to obtain the following rep- 
resentation for large m (Proposition 2.5) 

(^o,«"^(0)V'o) = (V-cPMi • • •M^_i_iiV^(0)PV'o) . (1.7) 

Here, P is the orthogonal projection onto Hs, along i/^c, projecting out the degrees of 
freedom of C. The Mk are effective operators which act on TCs only, encoding the effects of 
the interactions on the system S. They are called reduced dynamics operators (RDO), and 
have the form 

where is an (unbounded, non-normal) operator acting on Hs^'Hsi^, satisfying e'^^'^^e"'*^- 
gitLfc^g-itLfc fQj. Ae^s® ^£k, and X^Vs <^ V'S^ = 0-^ The operator N^{0) acts on 
Hs and has the expression (Proposition 2.4) 

Ar„(O)^0 = Pe^^—'-^— ' • • • e'^'^^^iAg 5^))e~'^"'-^'" • • • e-'^—'-^—'^o- (1-8) 

The asymptotics m ^ oo of (1.7) for identical matrices = M has been studied in [7]. 
In the present work we consider the to be random operators. We allow for randomness 

^These are the defining properties of Ki^; lias an explicit form expressible in terms of the modular 
data of (9Jt5 (gi ftHfj. , ips ® i'Sk), see Section 2.2. 
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through random interactions (interaction times, interaction operators) as well as random 
initial states of the £m (random temperatures, energy spectra, etc). 

Let {Q.,J^,^) be a probability space. To describe the stochastic dynamic process at 
hand, we introduce the standard probability measure dP on Oext '■= ^* , 

dP = nj>idpj, where dpj = dp, Vj G N*. (1.9) 

We make the following randomness assumptions: 

(Rl) The reduced dynamics operators are independent, identically distributed (iid) 
random operators. We write Mj. = M(a;jfc), where M : J7 — ^ B{C^) is an operator 
valued random variable. 

(R2) The operator N„i{0) is independent of the with l<k<m — I — 1, and it has the 
form N{oOm-i, ■ ■ ■ ,^m+r)-, where N : $7''+^+^ B{C'^) is an operator valued random 
variable. 

Since the operator describes the effect of the k-th interaction on S, assumption (Rl) 

means that we consider iid random repeated interactions. The random variable N in (R2) 

does not depend on the time step m. This is a condition on the observables, it means that 

(j) . 

the nature of the quantities measured at time m are the same. For instance, the in 
(1.6) can represent the energy of £m+j, or the part of the interaction energy Vm-\-j belonging 
to £m+ji etc. Both assumptions are verified in a wide variety of physical systems: we may 
take random interaction times = r(a;fc), random coupling operators Vfe = V{uk), random 
energy levels of the Ej. encoded in L^^ = Lsicok), random temperatures /Jg^^. = Psi^k) of the 
initial states of £k, and so on; see Sections 6 and 7 for concrete models. 

1.2 Main results 

Our main results are: the existence and identification of the limit of infinite products 
of random reduced dynamics operators; the proof of the approach of a random repeated 
interaction system to an asymptotic state, together with its identification; the analysis of the 
macroscopic properties of the asymptotic state; explicit expressions of that state for spin- 
spin and spin-fermion systems. We present here some main results and refer to subsequent 
sections for more information and for proofs. 

- Ergodic limit of infinite products of random operators. The asymptotics of the dy- 
namics (1.7), in the random case, is encoded in the product 

M{uJi) ■ ■ ■ M{Um-l-l)N{uJrn-l, ■ ■ ■,i^m+r)- 

It is not hard to see that the spectrum of the operators M(u;) is contained inside the closed 
complex unit disk, and that M{(jo)ips = i^S (see Lemma 2.3). 

Definition 1.1 Let M.(^e) denote the set of reduced dynamics operators whose spectrum on 
the complex unit circle consists only of a simple eigenvalue 1 . 

The following is our main result on convergence of products of random reduced dynamics 
operators (see also Theorem 3.3). We denote by E[M] the expectation of M(a;). 
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Theorem 1.2 (Ergodic limit of infinite operator product) 

Suppose that p(M(cl') G ^{e)) / 0- Then K[M] £ A^(e)- Moreover, there exists a set 
Q C ri^* of probability one s.t. for any uJ = (a;n)neN £ ^, 

1 " 

lim - y2M{u;i)---M{tOn)N{uJn+i,...,uJn+i+r+i) = \i^s){0\ E[N], 

i^— »0O 1/ ^— ' 

n=l 

where 6 = -Pi Ejjy^^jV'Sj *s i/ie (Riesz) spectral projection of X associated to the eigenvalue 
\, and * denotes the adjoint. 

- Asymptotic state of random repeated interaction systems. We use the result of The- 
orem 1.2 in (1.7), where we replace a"* by the random dynamics, denoted a^. It follows 
that the ergodic limit of (1.7) is ^+(E[iV]), where 

Q+{A):={e,A^s). A ems. (1.10) 

A density argument using the cyclicity of the reference state tpQ extends the argument 
leading to (1.7) to all normal initial states q on SDt. 

Theorem 1.3 (Asymptotic State) Suppose that p(M(a;) G M.{e)) 0- There exists a 
set Q, C r^^* of probability one s.t. for any cJ G for any instantaneous observable O, 
(1.6), and for any normal initial state q, we have 

lim iX]^(aL«(0))=^p+(E[iV]). (1.11) 

m=l 

- Macroscopic properties of the asymptotic state. Since we deal with open systems, 
it is generally not meaningful to speak about the total energy (which is typically infinite). 
However, variations (fluxes) in total energy are often well defined. Using an argument of [7] 
(see also [6] for a heuristic argument based on the hamiltonian approach) one shows that 
the formal expression for the total energy is constant during all time-intervals [Tm-iiTm)^ 
and that it undergoes a jump 

j{m,u) := a^{ViiOrr.+i) - VM) (1.12) 

at time step m. The variation of the total energy between the instants and m is then 
AE{m,uJ) = X^feLi j(^)^)- The relative entropy of g with respect to qq, two normal states 
on SDT, is denoted by Ent{Q\go). Our definition of relative entropy differs from that given 
in [5] by a sign, so that in our case, Ent(^|^)o) ^ 0. For a thermodynamic interpretation of 
entropy and its relation to energy, we assume for the next result that ipg is a {jSg, Q!^)-KMS 
state on dJts, and that the ip£^ are {Ps^, a£;^)~KMS state on Tls^, where Ps is the inverse 
temperature of S, and are random inverse temperatures of the Sm- Let be the state 
on 9Jt determined by the vector '4's ® '4'c = '4's ®m V'£m ■ The change of relative entropy is 
denoted AS'(m,cZJ) := Ent(^ o a'^\QQ) — Ent(^|^o)- 
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Theorem 1.4 (Energy and entropy productions, 2'^'* law of thermodynamics) Let 

Q be a normal state on Tl. Then 



lim J ^E(^^^) \ ^, dE+ = ^+fE[P(L5 + y-e-^(L5 + F)e--V]) 



a.s. 



lim M!!^ =. d5+ = ^+(E[/3£P(L5 + F-e'"^(L5 + ^)e-^"V]) cl.s. 

We call dE^ and dS+ the asymptotic energy- and entropy productions; they are independent 
of the initial state g. If is deterministic, i.e., ZJ -independent, then the system satisfies 
the second law of thermodynamics: dS-^- = PsdE^. 

- Explicit expressions for asymptotic states. We apply our general results to spin-spin 
and spin-fermion systems, presenting here a selection of results, and referring the reader to 
Sections 6 and 7 for additional results and more detail. 

Spin-spin systems. Both S and £ are two-level atoms with hamiltonians hs, hs having 
ground state energy zero, and excited energies Es and Eg, repectively. The hamiltonian 

describing the interaction of S with one £ is given by h = hs -\- hg -\- \v, where A is a 
coupling parameter, and v induces energy exchange processes, 

V := as ® a*g -\- a% ® as ■ (1-13) 

Here, a# denotes the annihilation operators and the creation operators of ^ = S,£. 
The Gibbs state at inverse temperature (3 is given by 



^^\^''*'^\ where = TV(e-^'^#). (1.14) 



We take the reference state to be = ^5 ®m>i '4^£m,0mi where ^5 is the tracial state on 
5, and ip£ra,i3£ is the Gibbs state of £m (represented by a single vector in an appropriate 
"GNS" Hilbert space, see Section 6). 

The following results deals with three situations: 1. The interaction time r is random. It 
is physically reasonable to assume that t{uj) varies within an interval of uncertainty, since it 
cannot be controlled exactly in experiments. 2. The excitation energy of £ is random. This 
situation occurs if various kinds of atoms are injected into the cavity, or if some impurity 
atoms enter it. 3. The temperature of the incoming atoms is random. This is physically 
reasonable since the incoming atom beam's temperature cannot be controlled exactly in 
experiments. 



Theorem 1.5 (Random spin-spin system) Set T :- 



2lT 



-^(Es-£;f)2+4A2 ■ 

1. Random interaction time. Suppose that Ps^ = (3 is constant, and that t{(jo) > is 
a random variable satisfying p (r ^ TN) 7^ 0. Then there exists a set Ct C VlP" of 
probability one, such that for all U ^ Q,, for all normal states g on Tl and for all 
observables A of S, 

lim - J2 ^«(^)) = Q/3',s{A), (1-15) 



jU— »0O a 

m=l 



with (3' = f3i:= f3Es/Es. 
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2. Random excitation energy of £. Suppose that r and = (3 are constant, and 
that Es{ijo) > is a random variable satisfying p (r ^ TN) / 0. (Here, T = T{oj) is 
random via E£{uj).) Then there exists a set C O^* of probability one s.t. for all 
aJ G for all normal initial states g onDJt and for all observables A of S, (1-15) holds 
with /?' = 132 := -Es^ log (2(1 - (1 - E[eo])-iE[(l - eo)(l - 2^-^^/^^ 5)] - l), 
and where 



eo 



[Es -Es- ^J{Es-EsY+A\^^ ^ + 4A2ei^V'(^s-S£)2+4A2 



(^Es -Es- ^{Es-EeY + AX^'^'^ + 4A2 



(1.16) 



3. Random temperature of £. Suppose that f3{(jj) is a random variable, and that r > 

satisfies r ^ TN. Then there exists a set Q, C 0,^* of probability one s.t. for alluJ G Q,, 
for all normal initial states g on and for all observables A of S, (1-15) holds with 

/?' = /?3 := -Es' log {nz-pi^)Es/Es,sr' - 1) • 

Remarks. 1. In the situation of point 1. of Theorem 1.5, we obtain the following sharper 
result than (1.15). There is are constant C, a > 0, and there is a random variable UQiu) 
satisfying E[e°"o] < 00 such that, for each aJ G \g{a'^{A)) - gp'^s{A)\ < Ce""", for ah 
n > no(aJ), all observables A and all normal initial states g. 

2. If Eg = Eg then /?i = (3. In the case of identical interactions (no randomness), the 
system 5 is therefore "thermalized" by the elements of the chain, a fact which was already 
noticed in [2]. One might expect that for a randomly fluctuating temperature of the f , 
the system S would be thermalized at asymptotic temperature equalling the average of the 
chain temperature. However, point 3. of the above theorem shows that this is not the 
case: the asymptotic temperature is in general not the average temperature. The random 
repeated interaction process induces a more complicated thermalization effect on S than 
simple temperature averaging. 

Spin-fermion systems. Let 5 be a spin- 1/2 system with Hilbert space of pure states C^, and 
Hamiltonian given by the Pauli matrix a^- We take the systems £ to be infinitely extended 
thermal fermi fields. They model dispersive environments. Let a{k) and a*{k) denote 
the usual fermionic creation and annihilation operators, and let a{f) = Jj^s f{k)a{k)d^k, 
^*(/) = /r3 /(^)<^*(^)d^fc, for square-integrable /. We take the state gp of £ to be 
the equilibrium state at inverse temperature f3. It is characterized by gp{a.* {f)a{f)) = 
(/, (1 + e^^y'f), where the h appearing in the scalar product is the Hamiltonian of a sin- 
gle fermion. We represent the one-body fermion space as f) = L2(M+, d//(r); g), where g is 
an auxiliary Hilbert space, and we take h to be the operator of multiplication by r G R"^.^ 
At each interaction step, S interacts with a fresh system £ for a duration r. The inter- 
action induces energy exchanges between the two interacting subsystems, it is represented 
by the operator XV, where A is a small coupling constant, and V = a'x<Si[a*{g) + a(g)]. Here, 
ax is the Pauli matrix and g = g{k) G L^(M^, d?k) is a form factor determining the relative 
strength of interaction between S and modes of the thermal field. We consider random 



* For instance, for usual non-relativistic, massive fermions, the single-particle Hilbert space is L^(M^, d'^fc) 
(Fourier space), and the Hamiltonian is the multiplication by \k\^ . This corresponds to g = L^{S^,dT,) 



(uniform measure on S^), and dy,{r) = ^r^^^dr. 
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interaction times of the form t{lo) = tq + o-(aj), where tq is a fixed value, and cr{uj) G [— e, e] 
is a random variable with small amplitude e. 



Theorem 1.6 (Random spin-fermion system) Assume that the form factor satisfies 
l2(k3 dSfc) < oo, and that p((7(u;) G |N — tq) / 1. There is a constant Aq > 
s.t. if < |A| < Ao, then Theorem 1.3 applies, and the asymptotic state q+, (1.10), has the 
following expansion: for any A G dJls, 



Q+{A) = q{a)AoQ + (1 - q{a))Aii + Ra,x{A), 



:i.i7) 



where A 



1-] 



{i,Aj), i,j = 0,1 and \0), |1) are the eigenvectors of az with eigenvalues 
±1. The remainder term satisfies \Ra^x{A)\ ^ C'||A||(e^ + A^), where C is independent of 
e, (T, A, A. 

The probabilities q(a) are given by 



Q!+ + a_ ' ■' Tq («+ + a_)^ 



where, with sinc(x) = sin(a;)/a;, 

|2 



d//(r) 



1 + e-/''- 



e '^'sinc^ 



{r T 2)ro 



+ smc 



(r ± 2)ro 



(1.18) 



^± = TO j d//(r) l^^^'^l [e-f^smc [(r ^ 2)ro] + sine [(r ± 2)ro]} 
V± = J Mr)j^^^ {e-^' cos [(r T 2)to] + cos [(r ± 2)to]} 



Expansion (1.17) shows in particular that to lowest order in A, the final state is diagonal 
in the energy basis. This is a sign of decoherence of S due to contact with the environment 
C. 



Organization of the paper. In Section 2 we cast the dynamical problem into a shape 
suitable for further analysis. Our main result there is Proposition 2.5. Section 3 contains 
the proof of Theorem 1.2, and in Sections 4 and 5 we present the proof of Theorems 1.3 and 
1.4, respectively. In sections 6 and 7 we present the setup and main results for spin-spin 
and spin-fermion systems. In particular, we give the proofs of Theorems 1.5 and 1.6. 



2 Repeated interactions and matrix products 

In this section, we link the repeated interaction dynamics to products of matrices. This 
reduction is a purely "algebraic" procedure and randomness plays no role here. Throughout 
the paper, we assume without further mentioning it, that 

(Al) dim = d < oo, and the reference vectors ip^ are cyclic and separating for Tl^ 
(# = 5 or Sm)- 
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Recall that cyclicity means that is dense in 7Y#, and separability means that 

= ^ = 0, G 9Jt#, and is equivalent to Tt'^ijj^ is dense in H^, where Tl'^ 
is the commutant von Neumann algebra of 

2.1 Splitting off the trivial dynamics 

We isolate the "free part" of the dynamics given in (1.2)-(1.4), i.e. that of the elements Sk 
which do not interact with »S at a given time step m. 



Proposition 2.1 For any m, we have 

U{m) = U- e-'rr.Lrn...^-ir,L, 

where 



Um = exp 



rn j—1 

j=i k=i 



and = exp 



j=i k>j 



(2.1) 
(2.2) 



are 



unitary operators which act trivially on Hs and satisfy U^ipc = V'C) Vm G N*. 



Proof. As the interaction Liouvillcan at time m, L^, and the free Liouvillean L^^^ commute 
provided k m, we can write successively 



-iT2L2 — g-iT-2i£i g-ir2L2 g-i'f"2 Z]fe>2 ■'^ffc 



(2.3) 



and then use this decomposition in (1.2). 



□ 



2.2 Choosing a suitable generator of dynamics 

We follow an idea developed recently in the study of open quantum systems far from equi- 
librium which allows to represent the dynamics in a suitable way [9, 7, 8, 12, 13, 14]. 
Let Jm and denote the modular conjugation and the modular operator of the pair 
{Tls (g) Tl£^,ips ipe^), respectively. For more detail see the above references as well as 
[5] for a textbook exposition. Throughout this paper, we assume the following condition on 
the interaction, without further mentioning it: 

(A2) A^/V^A"'/' ems® Ms^, Vm > 1. 



We present explicit formulae for the modular conjugation and the modular operator for 
the spin-fermion system in Section 7. The Liouville operator at time m associated to 
the reference state ^ps "X" ips^ is defined as Km = Ls + Ls^ + Vm 



'Jm'-^m ym'-^m <Jrr 



It satisfies ||e='=^-'^"' || < explHA^ Vm^m ||}- (In [9]; such operators are called C-Liouville 
operators.) The main dynamical features of Kjn are the relations 



Ae 



-itL„ 



JtK„ 



Ae"'*^-, G an^ (g) fmc,m >i,te 



(2.4) 
(2.5) 
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Relation (2.4) means that Km implements the same dynamics as L^- This is seen to hold by 

1/2 1/2 

noting that the difference Km — Lm = Jm^m Vm^m Jm commutes with all A £ (XiOJlc 
(since JdJlJ = 9H', as is known from the Tomita-Takesaki theory of von Neumann algebras, 
see e.g. [5]). The advantage of using K^, instead of is that e'*^"" leaves ips <X) ipe^ 
invariant. However, while is sclf-adjoint, K^ is not even normal and unbounded. 

We want to examine the large time behaviour of the evolution of a normal state g on 
9Jt, defined by o a™ (see (1-5)). Since a normal state is a convex combination of vector 
states, it is not hard to see that one has to examine the large time evolution of vector states 
only. More precisely, by diagonalizing the density matrix, we can write p = ^j>iPj\4'j){(l>j\j 
where pj > and ^j>iPj = 1, and where the are normalized vectors in Ti. If we can 
show that lim^^oo {4', a"^{A)4>) = ^'^(A) exists for any normalized vector (p & H, then any 
normal state satisfies 

hm e(a™(A)) = hm Y^pj (</.,-,«-( A) <^,) = Y^pjQ^.iA). (2.6) 

m— »oo m— »oo ' ■' ■' •' 

In other words, we only have to analyze vector states g{-) = {<p, -(p). If the asymptotic 
states do not depend on the vector (p, i.e. g^j, = g^, then any normal initial state g has 
asymptotic state by (2.6). The above argument works equally well if the pointwise limit 
m — > 00 is replaced by the ergodic limit. 

Next, since, by assumption (Al), ipo = ips'^ i'C, where ipc = ®m>i'^£mi cyclic for the 
commutant 9Jl' (which is equivalent to being separating for 93t), we can approximate any 
vector in Ti. arbitrarily well by vectors 

<P = ^Vo, (2.7) 

for some 

= 0^=1 i?; ®„>;v eOH', (2.8) 

with B'^ G 9Jt^, B'^ G 9Jt^^ (with vanishing error as ^ 00; see also [7]). Hence, we may 
restrict our attention to taking the limit m — >■ 00 of expressions 

(V^o, {B'Ya^{A)B'i;^) = (V^o, {B')* B' a"\A)^p^) . (2.9) 
2.3 Observables of the small system 

To present the essence of our arguments in an unencumbered way, we first consider the 
Heisenbcrg evolution of observables As G OTt^, and we treat more general observables in 
the next section. Consider expression (2.9). Using Proposition 2.1, we obtain 

a"^(^5®lc) = U{my {As®tc)U{m) (2.10) 
= (?7+)* e^^i-^i • • • e^^--^- {As ® Be) e"^^--^- • • • e-'^^^^ U+, 

where we made use of the fact that acts trivially on Hs- Due to the properties of the 
unitary i7+(m), specified in Proposition 2.1, and due to (2.4), (2.5), we have 

a'^iAs tc) ipo = {U^T • • • e'^--^"" {As ® Ic) e-'^-^- • • • e-'^'^'ipo 

= {U+y e^^i^i • • • e'^--^- {As ® lc)i^o- (2.11) 
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Let us introduce Pat = (g) l^-^ • • ■ Isj^ (g) Pip£^_^^ <^ ^^£n+2 ® ' ' ' ' where P^^^ = | t/)^- J (V'ffc | . 
From the definition of B' , (2.8), we see that (V'o|(-B')*-B' = {'^q\{B')* B' . Moreover, 
introducing the m-independent unitary operator 



:= exp 

we can write, for m > N, 



N-l N 
j=l k=j+l 



|Vo, {B')*B'{U+)*e'^^^^ ■ ■ ■ ei^^^^Pjve^^~+i-^^+^ • • • e^^^-^^C^s ]k;)V'o 
We define the projection 

p = ts<^mm, (2.12) 

and observe that 

Pjve^^^+i^^+i • • • e^'"'"-^'" (yl^ (g) ]lc)'0o = Pjve''"^+i^^+i • • • e'^'"'^'"P(^5 (g ]lc)'0o 

= Pe'^~+i^^+i • • • e'^'"-^-P(y45 (g llc)V'o- 

By a simple argument using the independence of the elements £k of C, we show exactly as 
in Proposition 4.1 of [7], that for any q > 1 and any distinct integers ni, - ■ ■ ,nq, 

p^iTniKn^ ^iT„2Kn2 . . . ^TnqK^^ p _ p^lTn^Kn-^ p^'Tn^K-n^ P . . . p^^'^nqKn^p (2.13) 

Therefore, introducing operators Mj acting on Hs by 

Pe^rjKj p = Mj®\il;c){il^cl or Mj ~ Pe'-'i^i P, (2.14) 
we have proven the following result. 

Proposition 2.2 Let As G Tis and (p = B'^q with B' as in (2.8). Then for any m > N 
we have 

{<P,a^{As®tc)<l^) (2.15) 

V^o, {B'YB'{U+r^^^''' ■ ■ ■ e'-^''^PMN+iMN+2 ■ ■ ■ M^{As ® lc)^o 

Proposition 2.2 shows how the large time dynamics of a repeated interaction system is 
described by products 

^rn = MiM2---Mm on Hs- (2.16) 

The main features of the matrices Mj, inherited from those of e^'^^^J, are given in the 
following lemma. 

Lemma 2.3 ([7], Proposition 2.1) Assuming (Al), we have Mjif^s = i^s, for all j G N*. 

Moreover, to any 4> G TLs there corresponds a unique A G OK^ such that (j) = Alps- •= 
II^IIbcHs) defines a norm, on 7is, and as operators on Hs endowed with this norm, the Mj 

are contractions for any j G N* . 

Remark: It follows from Lemma 2.3 that the spectrum of Mj lies in the closed complex unit 
disk, and that 1 is an eigenvalue of each Mj (with common eigenvector ips)- 
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2.4 Instantaneous observables 

So far, we have only considered observables of the system S. In this seetion, we extend the 
analysis to the more general class of instantaneous observables, defined in (1.6). Those are 
time-dependent observables, which, at time m, measure quantities of the system S and of a 
finite number of elements £k of the chain, namely the element interacting at the given time- 
step, plus the I preceding elements and the r following elements in the chain. Physically 
important instantaneous observables are those with indices j = —1,0: they appear naturally 
in the study of the energy exchange process between the system S and the chain (see Section 
5); they also appear in experiments where one makes a measurement on the element right 
after it has interacted with S (the atom which exits the cavity) in order to get indirect 
information on the state of the latter. 

The Heisenberg evolution of instantaneous observables is computed in a straightforward 
way, as for observables of the form As (8) Ic- We refrain from presenting all details of the 
derivation and present the main steps only. Let 

a™'"(5) := e'^^T=nrj)Ls,Se--^^ET=nrj)Le, ^ ^ < m, (2.17) 

denote the free evolution from time n — 1 to m of an observable B acting non trivially 
on TCci^ only, with the understanding that a^'" equals the identity for n > m. With this 
definition and (2.2), we get 

iU-nAs^^^^_iB(jJ)U- 

= ^5®,^=-,a-;7^-+\5C^)). 

Hence, 

a"^{As^^j=_iB(i^) (2.19) 
= (?7+)*e'^i-^i • • • ^^"^^"^{As a2fj^^^\B^J^^))e-'^^^^ ■ ■ ■ e-'^'^'U+. 

Consider a vector state {(p,-(p), where (p is given by (2.7). We proceed as in the previous 
section to obtain 

{4>,a"^iAs0'j=_i B^i))^) 
= (sVo, S'(C/+)*e^^i^i • • • e'^-^- (^5 ^'j=-i a';;,'++^+\B(JJ)) e"^^-^ 
= (sVo, S'(C/+)*P^e-i^^ • • • e--^- {As a:::;7^+'(5^)))^o) 
The vector to the right of (U^)* can be further expanded as 

giriKi . . . ^imK^p^^iTj,+^K^+, . . . ei^-^-(^5 (S^)))V'0 

= e^riK, . . . gir^if^vpe'^^+i-^^+i . . . ^^"^^^{As a^^^^^^{B^J,^))i;Q 
= e^^i-^1 • • • e^^^-^^PMjv+1 • • • Mm-i-i x 

^P^Tr^-iKr^-i . . . ei^-^-(^5 a';;^'^+^+\B^i^))i,o, (2.21) 




(2.20) 
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where P has been defined in (2.12), and where we have proceeded as in the derivation 
of (2.15) to arrive at the product of the matrices M^v+i • • ■ Mjn-i-i- We now define the 
operator = Nm{0), see (1.6), acting on Tig by 

(Nrni's) ® i'c := Pe--'^-' • • • e---^-(^5 ^ (i?^^i^))V'o. (2.22) 

We will also denote the l.h.s. simply by Nmtpo- The operator depends on the instanta- 
neous observable, Nm = Nm{As, Bm^\ ... , Bm). It can be expressed as follows. 
Proposition 2.4 Let a"*'" denote the dynamics from time n to time m, i.e., 

«"*."(.) = U{m,n)* ■ U{m,n), 
where U{m,n) = U{m)U(n)*, and U{m) is given in (2.1). Then we have 

iV^Vo = Pa^^^-'-\As(^]^_iB^i)i;o (2.23) 

k=i 

Proof. The second equality is clear, since the dynamics involves only the with indices 
k < m. To prove the first equality, we use the properties of the operators Kj and the 
definition (2.22) to see that 

iv^^o = Pe^^— '•••e'^'"^-(^5<8,'=-i«:::^;+'+'(^^^)) X (2.24) 

Next, we write the oi^^^^'^^ in terms of the generators Lg^^., see (2.17), 

Inserting this expression into (2.24) we can distribute the generators L£^_^. among the 
propagators in (2.24), and we see that 

Nmipo = Pe^'^—'^— ' • • • e^^'"^'"(yl5 (g)^=_; B^))e~'^'"^'" • • • e-^'^—'^—'^o, 

where the Lk, (1-2), give the full dynamics. □ 

Finally, can be defined on all of Hs in the following way. Prom Proposition 2.4, it 
is immediate that for all observables ^4^ in the commutant StJl^, we can set NmA'^ipo := 
A'^NmipQ. Since Tl'^ips = Hs (separability of ■^5), is defined on all of Hs- We have 
proven the following result. 

Proposition 2.5 Let O be an instantaneous observable, (1.6), and let cf) = B'i/jq with B' 
as in (2.8). Then we have for any m > N + I + 1 

(0, a'"(O)0) = (Vo, (S')*S'(c7+)*e-i^i • • • e-^^^PM^+i • • • M^_i_,Nm{0)^o) , 

where the Mj are defined in (2.14), and Nm{0) is given in (2.22). 

To understand the large time behaviour of instantaneous observables, we study the 

n ^ 00 asymptotics of products 

'^nNn+l+l= MiM2---MnNn+i+i ou Hs, (2.25) 

where A^^+;+i involves only quantities of the systems S and ffc, with k = n+l, . . . , n+i+r+1. 
The numbers l,r are determined by the instantaneous observable O (1.6). 
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3 Proof of Theorem 1.2 



According to Proposition 2.5, the large time dynamics is described by products of operators 
of the form (2.25), in the hmit n oo. We will use in this section our basic assumptions 
(Rl) and (R2), saying that the Mj form a set of iid random matrices, and that A^„+;+i is a 
random matrix independent of the Mj, j = 1, . . . , n. In this section, we review results of [8] 
on products of the form Mi • • • M„, and we extend them to products of random matrices of 
the form (2.25). Our main result here is Theorem 3.3. 

3.1 Decomposition of Random Reduced Dynamics Operators 

Let Pij denote the spectral projection of Mj for the eigenvalue 1 (c.f. Lemma 2.3) and 
define 

^r-=Pi,j^S, Pj:=m{ikj\, (3.1) 

where P^j is the adjoint operator of Pij. Note that {ijjj\tps) = 1 so that Pj is a projection 
and, moreover, Mji/;j = il^j. We introduce the following decomposition: 

Mj := Pj + QjMjQj, with Qj = l-Pj. (3.2) 

The following are basic properties of products of operators M^. 

Proposition 3.1 ([8]) We define Mq^ := QjMjQj. For any n, we have 

Mi---Mn = m{0n\ + Mq,--- Mq^, (3.3) 

where 

9n = ^n + M*Qjn-l + --- + M*Q^---M*Q,^r (3.4) 

= M*---M2>i (3.5) 

and where (V'5)^n) = 1- Moreover, there exists Cq such that 

1. For any j G N*, ||P,-|| = < Co and \\Qj\\ < 1 + Cq. 

2. sup{||Mq,^Mq._^ . . . Mq.JI, n G N*, ju G N*} < Co(l + Co). 

3. For any n G N*, ||6'„|| < C^. 

Typically, for matrices M^ G M.(^e) (recall Definition 1.1), we expect the first part in 
the decomposition (3.3) to be oscillatory and the second one to be decaying. 

3.2 The probabilistic setting 

We use the notation introduced at the end of Section 1.1. Let us define the shift T : fiext — 

^ext by 

{Tu;)j = LOj+i, V a; = (a;j)jgN e ^^ext- (3.6) 

T is an ergodic transformation of Jlgxt- The random reduced dynamics operators are char- 
acterized by a measurable map 

n3ui^ M{loi) G Md{C), (3.7) 
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where the target space is that of all d x d matrices with complex entries, d being the 
dimension of Tls- With a slight abuse of notation, we write sometimes M{iJ) instead of 
M{uJi). Hence, for any subset B C Md{C), p(M(uJ) G B) = p{M-^{B)) = ^^-"^{3) dp('^)> 
and similarly for other random variables. According to (Rl) the product (2.16) is *n(^^) '■= 
M{uJi)M{u2) ■ ■ ■ M{un) = M{T^uj)M{T^uj) ■ ■ ■ MiT'^-^uj). 

In the same way as in (3.1), we introduce the random variable ip{ijOi) G C defined as 

V'(aJ) := Pi(aJ)*^5, (3-8) 

where Pi(aJ) denotes the spectral projection of M(tJ) for the eigenvalue 1, and where * 
stands for the adjoint. We decompose 

M{uJ) := \iPs){i^m + Mq{lJ) = P{oj) + Mq{lJ) (3.9) 

as in (3.2). Note that iP{lJ) and Mq{lJ) define bona fide random variables: uJ i-^ Pii'^^) is 
measurable since U i-^ M(tJ) is [4]. In the next section, we will consider the process (see 
(3.4), (3.5)) 

9n{iJ) = M*{T''~^u)M\T'^-'^u)---M*{Tu)xp{u) (3.10) 

n 
3=1 

Note that 0„ is a Markov process, since 0„_|.i(a;) = M*(a;„+i)6'„(ZIJ). 

Finally, the operators A''^ = A^(0), given in (1.8) and Proposition 2.4, have the form 

Ar^(O) = iV(a;^_,, . . . ,ujm+r) = N{T^-^-^iJ), (3.11) 

see also condition (R2). 

3.3 Convergence results for random matrix products 

We have pointed out after Lemma 2.3 that the spectrum of any RDO lies inside the complex 
unit disk, and 1 is an eigenvalue (with the deterministic, i.e., clJ-independent, eigenvector 
^s). The following result on the product of an iid sequence of RDO's is the main result of 
[8]. 

Theorem 3.2 ([8]) Let M (to) be a random reduced dynamics operator. Suppose that p{AI{uj) G 
■^(E)) / 0- Then we have E[M] G A^(e). Moreover, there exist a set Q,i C O,^ with 
P(J7i) = 1, and constants C,a > 0, s.t. for any oJ G fii there is an no(w) so that 

||MQ(a;i) • • • MQ{un)\\ < Ce""", for all n > nQ{U), and (3.12) 
1 

lim -y^eJoJ) = 9. (3.13) 

n=l 

Also, no(w) is a random variable satisfying E[e""'°] < oo, and 

e = il- E[MQ]*r' = P*E[MiE[^] = Pi,EiM]^s. (3.14) 
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As a consequence, 

1 " 

lim - y2MiiVi)---M{iVn) = \i^s){0\=Pi,ElM]- (3.15) 

I/— >0O 1/ ' ' ' 

n=l 

Remark. In the setting of Theorem 3.2, if not only M{u>), but also M*{uj) has a deter- 
ministic eigenvector with eigenvalue 1 (denoted and normalized as (ip'^jips) = 1), then 
9 = il^Q and one can sharpen (3.15) as follows (see Proposition 3.1 and equation (3.12)): 
There are constants C, a > 0, and there is a random variable rao(^) with E[e""''] < oo, s.t. 
for all aJ G Qi and all n > no(aJ), we have ||M(a;i) • • ■ M{(jjn) - \'4^s){il^*s\ II ^ C'e~°^". 

While this result allows us to study the large time behaviour of observables of the 
small system 5 (see Section 2.3), in order to study the physically relevant instantaneous 
observables, we need to understand products of the form (2.25). In our probabilistic setting, 
they read M{uj{) ■ ■ ■ M{LOn)N{T"'Lj). 

Theorem 3.3 Let M{lo) be a random reduced dynamics operator and let N{lJ) be a random 
matrix, uniformly bounded in u. Suppose that p{M{uj) G A^(£;)) ^ 0. Then there exists a 
set ^2 C Q^* s.t. ^(^2) = 1 and s.t. for any ZJ E ^2, 

1 " 

lim -V M(a;i)---M(cc;^)iV(r"uJ) = \i^s){0\ ^N]. (3.16) 

n— >oo 1/ ' 
n=l 

Remark. In our dynamical process, N(lJ) depends only on finitely many variables 
Ujn-i, . . . jUJm+r, scc (3.11), SO measurability and boundedness of the random matrix N 

are easily established in concrete applications. 

Proof of Theorem 3.3. Using the decomposition (3.3) together with (3.12), it suffices to 
show that 

lim - y 7V*(r"cZ;)^„(tJ) = E[N]*0. (3.17) 

n=l 

We follow the strategy of [8] used to prove (3.13) of the present paper. Prom (3.4) we get 

ly V n—\ 

^ iV*(r"a;)e„(a;) = ^ ^ iV*(r'^cJ)M^(r'^-^a;) • • • M^(r^+ia;)V'(T^u;) 

n=l n=l j=0 

V p—k 

= ^^Ar*(r'=+^cj)M^(r'=+J-^cj)---M^(rJ+^cj)V'(T^cj). (3.18) 

k=l j=0 

Let us introduce the random vectors 

9^''\lJ) = N*{T^u)M^{T^-^u)M^{T^-'^u) ■ ■ ■ M^{'nj)'iP{T^iJ), (3.19) 
so that, by (3.18), 

u—k 



n=l k=l j=0 

oo p—k ^ oo 

= i;X{fe<.}E^^'^(^'^)-=^E5(^''^'^)- (3-20) 



fe=l j=0 k=l 
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For each fixed k, by ergodicity, there exists a set fi^j.) C Jl^* of probabihty one, such that, 
for all lJ G i^(fc), the following limit exists 

lim gik,v,Uj) = lim ] Y e^''\T^Uj f ~ ^ ^ 

i^— »oo i^— »oo U — k + \ V 

1 ^ 

= lim -— - y 9^^\T^U) = E[e^% 

j=0 

Therefore, on the set J7oo ■= ^keN^{k) oi probability one, for any k G N, we have by 
independence of the M{uj), 1 < j < k, and of N*{T'^uj), 

lim g{k,u,Lj) = E[0(*=)] =E[N*]E[M^f-^E[iP]. (3.21) 

It follows from Proposition 3.1, Theorem 3.2 and the boundedness of N{u;) that for Zo G 
^2 = ^inn^, we have \\9^''\T^i:j)\\ < Ce'*no(T%)g-a(fc-i)_ Therefore, for all u large enough, 
and for all 1 < k < u, 

^ v—k 

\\g{k,v,Lo)\\ < C-^e""°(^'^)e-"('=-^) < 2(:7E[e""o]e-"(*^-^\ (3.22) 

i=o 

where we have used ergodicity in the last estimate. Of course, the same upper bound (3.22) 
holds for k > u, since then g{k,u,uj) = 0. The r.h.s. of (3.22) is summable w.r.t. A; G N, 

so we can use the Lebesgue Dominated Convergence Theorem in (3.20) to conclude that, 
almost surely on ^2, hm^^oo ^ En=i iV*(T"cJ)^n(cJ) = E[N*] J2T=oHMq]'' E[iI;]. Relation 
(3.17), and thus the proof of the theorem, now follow from (3.14). □ 



4 Proof of Theorem 1.3 

Let be a normalized vector in Ti.. Fix e > and uJ G ^^ext- There exists a B' = B'{e,Io) G 
9Jl' of the form (2.8) (with A'' depending on e,uJ), s.t. 

||</)-SVo|| <e. (4.1) 

Here, both ^ and may depend on W. It follows that 

I (<^,a^(0).^) - <5Vo,<(0)SVo> I < 2e||0||. (4.2) 

Using Proposition 2.2 and Theorem 3.3, and that B' commutes with a^(0), we arrive at 
the relations 

lim lf](^o,(S')*5'<(0)Vo>= lim - (^o, (S')*5'a^(0)Vo> 



= lim - y (V'o, (B')*B'(c7+)*e^^('^i)'^('^i) • • • e^^('^~)^('^^)px 

^ m=N+l 

xM{uJN+l) ■ ■ ■ M(Um-l-\)N{Um-htOm-l+U ■ ■ . ,UJm+r)tpO^ 

= (V'o,(S')*S>o> {e,E[N{0)]i;s), (4.3) 
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for all tJ in a set Q2 of measure one. It follows from (4.1) that (1 — e)^ < {tpo, [B')* B'iJjq) = 
ll-^VolP < (1 + e)^- Since e is arbitrary, using the latter bound in (4.3) and taking into 
account (4.2), we conclude that (1.11) holds for any vector initial state q{-) = (0, -0). 
Finally, the argument leading to (2.6) shows that (1.11) holds for all normal initial states. 
The proof of Theorem 1.3 is complete. □ 



5 Proof of Theorem 1.4 

An easy application of Theorem 1.3 shows that for any normal initial state q, 

hrn^ = Q+ij+), (5-1) 

where 

j+ = E [PVP - Pe'^^Ve-'^^P] = E [P{V - a^{V))P] . (5.2) 

The energy grows linearly in time almost surely, at the rate dE+.^ In order to show the 
expression for d£'+ given in Theorem 1.4, it suffices to prove that Q+\E{P{Ls—a'^{Ls))P)] = 
0. 

Let ^ be a normal state. Although ^ dJl, still a''{Ls) — a^~^ (Ls) is an instantaneous 
observable belonging to 9Jt. This follows from e^'^kLk^gQ-^TkLk _ g ^ msf., which in 
turn is proven by noting that 

^iTkLki^Q-iTkLk _Ls= / ' e"-^'=[iLfc,L5]e-'*-^'=di = / ' e'*-^'=[iVfe,L5]e"'*^Mt, 

Jo Jo 

where [iVk,Ls] = -^e^^^Ffce""^^ |t=o eMs^ me,. 

As a consequence, a^{Ls) — {Lg) G and we can apply Theorem 1.3 to obtain 
lim^^oo ^ T7k=i Q {<^\Ls) - a^+\Ls)) = Q+{nP{Ls - a-{Ls))P\) a.s. On the other 
hand, we have that Y.'k=i Q [oi^i^s] - a^+^{Ls)) = {a^iLs) - a™+^(L5)), which tends 
to zero as m — >^ GO. This proves the formula for dE^ given in Theorem 1.3. 

Next we show the expression for dS-^- in Theorem 1.3. The following result is determin- 
istic, we consider Zo fixed and do not display it. 

Proposition 5.1 Let g be a normal state on m. Then for any m> 1, we have 
AS{m) := Ent{g o a'^lgo) - Ent{g\go) 

(m 
T^PSk (jik) + oc'-\Ls + 14) - a\Ls + Vk+i)) + Psia'^iLs) - Ls) 
k=l 

where the energy jump j{k) has been defined in (1.12). 

®The definition of AE+ differs from the one of [7] by a factor i: here d-E+ represent the asymptotic 
average energy production per interaction and not per unit of time. One could also study the average energy 
production per unit of time. It is easy to see that 

g(Ag(m,a;)) ^ dE+^ ^ ^ 
m— oo r(a;i) H h r(a;m) E[r] ' 
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Proof. The proof is similar to that one of Proposition 2.6. in [7]. Using the entropy 
production formula [10], we have 

AS{m) = e(a"'(^(3sLs + Y,f3£,Le,'^ - ^/J^.L^, - /^^L^j • (5.3) 

Clearly, the sums in the argument of g in the right hand side only extend from = 1 to 
k = m. We examine the difference of the two terms with index k. 

= (3£,a''{Lk)-P£,L£^-P£,a''{Ls + Vk) 

= P£,a^-\Lk) - P£,L£^ + P£j{k) - f3£,a\Ls + F^+i) 

= l3£,a^-\Ls + Vk) + f3£j{k) - ^£,a\Ls + T4+i), 

where we use a^{L£^_) = a^{L£^_) in the first step, a^{Lk) = a^~^(Lk) and (1.12) in the 
third step, and in the last one a'^''^{L£^) = L£^. □ 

By Proposition 5.1, we have for all m > 1, a; G figxt 
A5(m,a;) 1 ^ 1 



fc=i fc=i 

{a^-\Ls)-a\Ls)))+-Q{[3s{or{Ls)-Ls)). 



m — ' V V J J m 

k=l 



Using Theorem 1.3 we see that with probability one (and where M denotes the reduced 
dynamics operator) 

lim ^^i^^"^) = (eIPsM] E[PVP] - E[f3£Pa^{V)P]) + q+ (e[P£PVP] 

m— »oo m \ /V 

-E[P£M] E[PVP]^ + g+(E[p£P{Ls - a^Ls))P]) 
= g+(E[(3£P{{Ls + V)-a^Ls + V))P]y 

This completes the proof of Theorem 1.4. □ 



6 Spin-spin models and proof of Theorem 1.5 

In this section, we consider both S and £ to be two-level systems, with interaction given by 
(1.13). This is a particular case of the third example in [7]. The main results of this section 
have been anounced in [8]. 

The observable algebra for S and for f is 2I5 = Slg = M2(C). Let Es,E£ > Ohe the 
"excited" energy level of S and of £, respectively. Accordingly, the Hamiltonians are given 

by 


E£ ■ 



hs = 




Es 



and h£ = 
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The dynamics are given by Q'^g{A) = p;^t^s j^^-tths ^^^^ = c**'^^ Ae"**'*^ . We choose 

(for computational convenience) the reference state of £ to be the Gibbs state at inverse 
temperature /3, see (1.14), and we choose the reference state for S to be the tracial state, 
^o,5(^) = The interaction operator is defined by \v, where A is a couphng constant, 

and V is given in (1.13). The creation and annihilation operators are represented by the 
matrices 



" 


1 ■ 




" 


" 


and = 








1 






The Heisenberg dynamics of S coupled to one element £ is given by the *-automorphism 
group t ^ eith^Ae-^^^^ , A e ®^Sj hx = hs + hg + Xv. 

To find a Hilbert space description of the system, one performs the GNS construction 
of (2l5,£»o,5) and {^e,Q/3,s)i see e.g. [5, 7]. In this representation, the Hilbert spaces are 
given by Us = Ti-e = ® C^, the Von Neumann algebra by Ms = Ms = M2(C) C 
B{C'^ C^), and the vectors representing ^o,5 and Qjj^s are i^s = ^ (|0) (E) |0) + |1) (gi |1)) 

and # = (|0) ^ |0) + e-^-^g/^|l) |1)), respectively, i.e., wc have ^/3#,#(^) = 

{A® ]1)V'#), # = S,£, Ps = P, Ps = 0, and where |0) (resp. |1)) denote the ground 
(resp. excited) state of hs and hg. Finally, the Liouvillean L is given by 

L = (/l5 (g) ]lc2 - ]lc2 ^5) ® (1C2 ® %) + (]1C2 1C2) <^ (^£ ® - ® ^f) 

+X{as (8) ]lc2) <g {a} (g) lea) + X{as (g) V) (g {as (g Ic^). 

6.1 Spectral analysis of the reduced dynamics operator M 

The RDO M is defined by (2.14). However, in this example, where the hamiltonian h\ is 
explicitly diagonalizablc, wc shall use another expression for it, which may look less simple 
but has the advantage that it only makes use of the self-adjoint hamiltonian. Since is 
cyclic for SDt^ and Hs has finite dimension, G Hs, 3!^^ = A (g 1^2 G VJls such that 
(f) = Asil^s- It is then easy to see that 

M{A (g lc2)ips = {M{A) (g lc^)ips, (6.1) 
and where the map M acts on ^s and is defined as 

M{A) := Its (e'^^^ A ® 11 e"'^''^) , (6.2) 

where Trs{As (g As) := Qf3,s{^£)As denotes the partial trace over £. 

Similarly, if M* denotes the map dual to M, i.e. Vp,^ G M2(C), Tt{pM{A)) = 
Tr(A^*(p)^), then we have, for any density matrix p 

i{M*{p)r(E)l)^s = M*ip*(E)l)iJs- (6.3) 

In particular, the spectrum of the map Ai* is in one-to-one corrcspondance with the spec- 
trum of the operator M* (via complex conjugation), and if p is an eigenvector of Ai* for 
the eigenvalue 1 (which we know to exist), then the "corresponding eigenvector" of M* is 
V'5 = {p* ® T^)'^s- A simple computation shows that the four eigenvalues of hx are £'0+ = 0, 
Eq^ = Es + Es and 

Ei± = \{Es+ Es) ± ^v'(£^5-^£)2 + 4A2. (6.4) 
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The corresponding normalized eigenvectors are given by ■i/'o+ = |0) <8) |0), ipQ- = |1) (8) |1), 
and ipi± = ai±|l) (8) |0) + bi±\0) (8) respectively, where 

ai+ = I, oi+ = — I. (6.5) 

v/A2 + (E5-^i±)2 ^\^ + {Es-E^^Y 

We finally denote ao+ = 6o_ = 1 and oq- = 6o+ = 0- Inserting the spectral decomposition 
of h\ into (6.2) gives the following result. 

Lemma 6.1 For any A G 21, 



n,(T,n ,(T' 

X 



{a.nc0.n'c'\f^'\"l^\ + e ^^^hnchn'c'\^ " " '^'l) > (6-6) 



where n,n' G {0,1} anc/ a, a' G {— >+} a^c^ -^/J.g = Tr(e ^'^^). Similarly, for any density 
matrix p, 

M*{p) = Z-^l^ e^^^^--^"'-^') {an^-ar,^a'{n\pn) + e-^^^6na6n'a'(l - "V(l - "))) 

n,<T,n',(T' 

X (an(^a„v|n')(n| + 6n(^6„'cr'|l - - '^■1) • (6-7) 
The above lemma allows us to make a complete spectral analysis of M. 
Proposition 6.2 

1. The eigenvalues of M are l,eo,e_,e+ where cq is given in (1.16), 

[Es -Es- y/{Es - EeY + 4A2) ^ + /^x^e'^^/{Es~Ee?+i\'' 



e_ 



-Es- ^{Es-E£y + 4X^y + 4A2 



Moreover, the eigenstates of M* for the eigenvalues l,eo,e_,e+ are respectively = 
[e-^'^'s lc2)i;s where /?' := pEg/Es, (^o = |0) (8 |0) - |1) (g) = |0) |1) and 

</>+ = |l)®|0). 

^. The functions |eo(T)|, |e+(r)| and |e_(r)| are continuous and periodic of period T := 
= . Moreover, they have modulus strictly less than 1 if and only ifr ^ TN. 
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-^(£;5-E£:)2+4A' 

Remark 6.3 Since cq is positive, point 2. proves that 1 is a non degenerate eigenvalue for 
M.* if and only if r ^ TN, i.e. for all but a discrete set of interaction times. This condition 
agrees with the corresponding assumption of [7] in the perturbative regime. 

Proof. Point 2. follows from point 1. Point 1. is proven by direct computation using 
(6.4)-(6.5)-(6.7). □ 
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6.2 Proof of Theorem 1.5 



Point 2. of Proposition 6.2 shows that M G M(e) if ^-nd only if r ^ TN. Hence, for this 
spin-spin model, Theorem 1.3 applies if and only if p(t ^ TN) ^ 0, which is precisely the 
assumption we have in each of the three situations of Theorem 1.5. It remains to compute 
the asymptotic state g+ in each of these three situation. Using the complete spectral 
decomposition of M*{uj) (see Proposition 6.2), we compute explicitly its expectation E[M*] 
and then the spectral projection Pi^e[m*]- After computation, we get: 

1. Random interaction time: Pi,e[m*] = j^^^3i3/hg^ IV'g) (e~^'^-^ lea 'tps\, 

2. Random excitation energy of £: -Pi,e[m*] = \i's){PE ® Ic^ '^s\, where 



PE 



(1 - E(eo))-^E((l - eo)(l - 2Z,,]^))\ |0)(0| 



+ 



1 + (1 - E(eo))-^E((l - eo)(l - 2Zji ))! |1)(1|, 



3. Random temperature of £: Pi^e[M* 



V'^l) where P 



Combining these formulas with (1.10) give the various expressions for the asymptotic state 
Finally, when the interaction time r is random, the map M*(lu) has a deterministic 
eigenvector for the eigenvalue 1. This allows for stronger convergence results as mentioned 
in the Remark after Theorem 3.2. 



7 Spin-Fermion models and proof of Theorem 1.6 

We combine our convergence results with a rigorous perturbation theory in the coupling 
strength between S and C. We take <S to be a 2-level atom and the £ are large quantum 

systems, each one modeled by an infinitely extended gas of free thermal fermions. The ran- 
dom parameters are the temperature of the system S^, = Pj^^ , as well as the interaction 
time Tfe. 

The state space and the reference vector of S are 

Hs = C^(^C^ V5 = ^(|0)®|0) + |1)®|1)), (7.1) 

where {|0) = [1,0]"^, |1) = [0,1]'^} is the canonical basis of C^. (7.1) gives the GNS 
representation of the trace state on the algebra of complex matrices M2(C), ^Tr{As) = 
{il^s, {As (g) ls)ips), for all As G M2(C). The von Neumann algebra of observables repre- 
sented on Tis is thus Tls = M2(C)0l C ^(C^^C^). The Heisenberg dynamics of 5 IS given 
by e^*'^^ A^e^'*'^"'. The Pauli matrices cj^ and ax (the latter plays a role in the interaction) 
are 





" 1 " 




' 


1 " 




-1 




1 






(7.2) 
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On the algebra Tls, the dynamics is implemented as Tg{As (8) 1) = e^^^^{As ll)e 
with standard Liouville operator 

Ls = crz<^l-l<^az. (7.3) 

Note that Lgips = 0, as required in (1.1). It is easily verified that the modular operator As 
and the modular conjugation Js are given by 

A5 = 1(8)1, J5('0«'X) = (7.4) 

for vectors ■0) X ^ C^, and where the bar means taking complex conjugation of coordinates 
in the canonical basis. 

We now describe a single element £ of the chain, a free Fermi gas at inverse temperature (3 
in the thermodynamic limit. We refer the reader to [5] for a detailed presentation. Let f) and 
h be the Hilbert space and the Hamiltonian for a single fermion, respectively. We represent 
^ as f) = L^{M.'^,dfi{r);g), where g is an auxiliary Hilbert space, and we take h to be the 
operator of multiplication by r G R"*". (See also footnote 4 at the end of Section 1). The 
fermionic annihilation and creation operators a{f) and a* (/) act on the fermionic Fock space 
r_([)). They satisfy the canonical anti-commutation relations (CAR). As a consequence of 
the CAR, the operators a{f) and a*{f) are bounded and satisfy ||a*(/)|| = ||/|| where a* 
stands for either a or a*. The algebra of observables of a free Fermi gas is the C*-algebra of 
operators 21 generated by {a*(/)|/ G f)}- The dynamics is given by r£(a*(/)) = a*(c^*'*/), 
where h is the Hamiltonian of a single particle, acting on f). It is well known (see e.g. [5]) 
that for any P > 0, there is a unique (rp,/?)— KMS state gp on 21, determined by the two 
point function Qi3{a*{f)a{f)) = (/, (1 + e^^^)^^f). Let us denote by Jlf the Fock vacuum 
vector, and by N the number operator of r_(()). We fix a complex conjugation (anti-unitary 
involution) / ^ / on which commutes with the energy operator h. It naturally extends 
to a complex conjugation on the Fock space r_([)) and we denote it by the same symbol, 
i.e. $ ^ 

The GNS representation of the algebra 21 associated to the KMS-state is the triple 
(W£,7r/3,V'£) [1] where 

= r_(l))®r_(f)), i>s = ^f®^U (7.5) 

and 



Ma\f)) = a* (^/) « 1 + (-1)^ ^ a (^/) =: a^(/). 



(7.6) 



The von Neumann algebra of observables for an element £ of the chain is dJls = 7r/3(2t)", 
acting on the Hilbert space He- The dynamics on TTp{^2i) is given by Tl{7rf){A)) = Trp(T^{A)), 
it extends to Wis in a unique way. The standard Liouville operator is given by 

Le = dr{h) 1 - 1 (8 dr(/i). (7.7) 

Note that Lg^s = 0. Finally, the modular conjugation and the modular operator associated 
to {M£,il;£) are 

j£($(8*) = (-l)^(^-i)/2*(8(-l)^(^-i)/2$, A£ = e-^-^^. (7.8) 
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The combined, uncoupled system has product structure, with Hilbcrt space Tis ® Tie, 
algebra Tts dJls , reference state tps ® ipe- The uncoupled dynamics is generated by the 
Liouville operator 

Lo = Ls + Le. (7.9) 

We now specify the interaction between the small system and the elements of the chain. 
Let g e i) he a form factor. The interaction operator is given by 

V:=a^0 ]lc2 ^ Mg) + a}{g)) (7.10) 

(where ax is defined in (7.2)). It produces energy exchange processes between S and £. 
Using (7.4), (7.8), one readily calculates 



(A^ ® A£)VV(A5 A^)-i/2 
= cr j; (g) lc2 $ 



1 \ ^ f 

g] ®l + a[ 



+(-1)^ ® a* g + (-if ® a g 



(7.11) 



We assume that e^'^^^g G I). Then (7.11) shows that {As (8) As^/'^ViAs (8) Ag)-^/^ ^ 
<8> 9?l£, i.e., Condition (A2) of Section 2 is satisfied. 

Theorem 7.1 (Convergence to asymptotic state) Let < Tmin < ^max < oo and < 
Praax < OO be given. Let r : i-^ [TmimTmax] O'nd /? : $7 1-^ (0, /3max] be random variables. 
Suppose that ||(1 + e^™'"'''/^)^!! < oo, and that there is a 6 > such that 

p(dist(r,fN) >(5) 7^0. (7.12) 

Then there S, and on the form factor g, 

s.t. if < |A| < Ao, then p(M{uj) € Mi^e)) > 0- -^'^ particular, the results of Theorem 1.3, 
applied to the spin-fermion system, hold: the system approaches the repeated interaction 
asymptotic state q+, defined in (1.10). 

Proof. We expand the operator M in a power (Dyson) series in A: 

M = e'^'^^P (7.13) 
+ V A^" / dti • • • / dt2„ e''^^-5Pe"2"-^ot^e"'*2n-ffo . . . ^ihKo^^^-itiKo 

n>l 



where only the even powers appear since the interaction is linear in creation and annihilation 
operators, and P projects onto the vacuum. W is the operator 

W = V- JA^/VA-^/V, (7.14) 

where V is given in (7.10) (see also (7.11)), and J = ® Js, A = As (g) A^- are the modular 
conjugation and the modular operator associated to {Tts ® Tl£,ijjs ® ips), see also (7.4), 
(7.8). Using the Canonical Anticommutation Relations, one easily sees that ||a*(5)|| = H^H 
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(independent of /3; see (7.6) for the definition of the thermal creation and annihilation 
operators). Using (7.14) and (7.11), it is easy to find the upper bound 

IIT^II < 3|| (1 + 6^^^/2)^11 _ ^7J5) 

We apply standard analytic perturbation theory to the operator (7.13). For A = 0, the 
eigenvalues of M, {Ije^^'"^}, lie apart by the distance 

ro(r) :=min{2|sin(2r)|,2|sin(r)|}. (7.16) 

(Note the spectrum of Ls is {—2,0,0,2}, c.f. (7.3).) We assume that the interaction time 
is such that ro(r) is strictly positive. Below, this condition appears as dist(r, |N) > S > 0. 

The following result gives an estimate of the eigenvalues of M, which will be needed in 
verifying that M is in the family M(^e)- 

Proposition 7.2 Suppose that |A| < \ro{T). Denote by l,eo,e± the four eigenvalues of 
M. We have 

eo = 1 - A^T^a + £0, (7.17) 



e± = e±2- 



2^2 



l-^a±iAV /dM(r)||5(r)f,x 



/ Mr)\\g{r 

-L r\\ \ 1 

+ £±, (7.18) 



/ 1 — sinc(r(2 — r)) 1 — sinc(T(2 + r)) 
2^ + 2+7 

where sinc(a;) = sin(x)/x and where 

a = I dMr)||5(r)||^(sinc2[^^] +sinc2[^^]). (7.19) 

The error terms # = 0, ±, satisfy the bound 

\s#\ < 12X^r^\\Wfcosl.H\X\rm)[l+'-^^ (7.20) 

Proof. Expansions (7.17), (7.18) of the eigenvalues have already been calculated in [7], Sec- 
tion 4.8, but the error estimate (7.20), allowing the control of r, /3, has not been given there. 
This error estimate is obtained by performing perturbation theory in a straightforward, but 
careful fashion. One proceeds as in [11], Chapter II.2. □ 

By knowing this expansion of the eigenvalues of M, we can impose a smallness condition 
on A which guarantees that the eigenvalues e# have modulus strictly less than one, which 
is equivalent to saying that M £ Mi^e)- 

Proposition 7.3 Suppose that Tmin < r < Tmax; /? < Pmax, o-nd that dist(T, ^N) > 5, 
for some constants < Tmm < ^"max cind /?max; (5 > 0. Then there is a constant Aq > 0, 
depending on Tmin, Tmax, /Smax, o,s well as on the form factor g , s.t. «/0 < |A| < Aq, then 



\e#\<l-^a<l, (7.21) 



# = 0, ±. In particular, M G M.(e)- 
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End of proof of Theorem 7.1, given Proposition 7.3. Fix r^mn-, 7"max5 /^max and (5, and 

suppose that (7.12) holds. Denote by Q,' the set of lo for which dist(r(u;), ^N) > 5. Then 
p(ri') / 0, and for each uj e O', we have M{u)) G by Proposition 7.3. Consequently, 

p(M(a;) G > p(J^') > 0. □ 

Proof of Proposition 7.3. We impose conditions s.t. the three eigenvalues given in 
(7.17), (7.18) have modulus strictly less than one. We have 

|eo| < 1 - -^a, (7.22) 

provided 

kol < (7.23) 

Next, since e_ is the complex conjugate of e+, it suffices to consider the latter. We write, 
with obvious identifications in (7.18), e+ = e^''^[l — x + \y + £+]. We have 



|e+| <\l-x + \y\ + \e+\ = yjl - 2x + x"^ + y"^ + \e+\. (7.24) 



Since + is just the square of the modulus of the second order (A^) contribution to the 

1 ' 
ro(T) - 



eigenvalue e+, it is easy to see that x^ + y"^ < A^t^||W||^(1 + 7777)^- We now impose the 



condition 



1 ^2 



2AV||VF||M 1 + — - <a, (7.25) 
^o(r)/ 



which implies that x + y < x. Combining this latter inequality with (7.24) gives 



|e+| < v^T^+|£+| < l-a:/2 + |£+|. (7.26) 
Finally we impose the condition 

\ 2 2 

|£+|<x/4 = -^a, (7.27) 

so that we get from (7.26) 

\2 2 

|e+|<l-x/4 = l ^a. (7.28) 

8 

This last bound, combined with (7.22), proves that (7.21) holds, provided the conditions 
(7.27), (7.25) and (7.23) are imposed. Taking into account the bound (7.20), we see that a 
sufficient condition for (7.27), (7.25) and (7.23) to hold is that 

mX\-\\Wtcos^.\\X\r\m)\l+'-±^^^ (7.29) 

L ro(r) J 

One may now use (7.15), (7.16), to find a constant Aq, depending only on the parameters as 
stated in the proposition, s.t. if jAj < Aq, then (7.29) holds. (Note that a, (7.19), does not 
depend on /3, and the minimum of a, taken over r > varying in any compact set, must be 
strictly positive.) This completes the proof of Proposition 7.3, and hence that of Theorem 
7.1. □ 
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7.1 Proof of Theorem 1.6 



Since E[M] G ■M.(e) (by Theorems 7.1 and 1.2), 1 is a simple eigenvalue of E[M]. Let ip"^ 
denote the unique vector invariant under E[M*], normalized as ('i/'J,^^) = 1, where ips is 
given in (7.1). We have Pi,e[m] = \'>Ps){tPs\^ ^"^^ t'^^^' (3-14), 9 = ip"^. To calculate ip^, 
we note first that for any u>, M{oj) is block-diagonal: 

Lemma 7.4 Let Pq = |0)(0| + |1)(1| be the spectral projection of Ls associated to {1}. 
The operator M {oj) leaves the subspace RanPo invariant. In the ordered orthonormal basis 
{|0),|1)} o/RanPo; we have the representation 



PoM{u)Po = 1- AV2(w 



— Q!+(a;) a+ii-^) 



+ 0(A^ 



(7.30) 



where the a±(w) are given by (1-18) with t{ijj) replaced by tq. The remainder term is 
uniform in r varying in compact sets. 

Proof of Lemma 7.4. As explained at the beginning of the proof of Theorem 7.1, only 
even powers of the interaction are present in the Dyson series expansion for M, (7.13). 
It follows from (7.10) and (7.14) that each term in the Dyson series (7.13) leaves RanPo 
invariant; this is so because the operator ax shows up an even number of times, and cr^ |0) = 
|1) and (Ta;|l) = |0). The calculation of the explicit form (7.30) is not hard. This concludes 
the proof of Lemma 7.4 □ 

The expansion for M* and hence of E[M*] in powers of A follow directly from (7.30). One 
then performs an expansion in powers of a and finds for the 0(A^)-term: 



a- 

— O- 



- 2X^E[a] 







-C- 


e+ 


V- 


-v+ 


-V- 


v+ 



The following expansion of the invariant vector ip^ follows: 



1 
72 



+2E[a 



-|0) |0) + 



-|1)®|1) 



2....^(|0>«|0) 



+ 



T§{a+ + a-) 



+4(EH)2(e_+e4 



r^(a+ + a 



+E[a' 



ro(Q;+ + a_)2 



3(|0)®|0)-|1)®|1)) 
(|0)® |0) -|1)® |1)) +0(a3) + 0(A2) 



(7.31) 



Formula (1.17) now follows directly from (7.31) and (1.10). This concludes the proof of 
Theorem 1.6. □ 
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